Transmission through a biased graphene bilayer barrier 
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We study the electronic transmission through a graphene bilayer in the presence of an applied bias 
between layers. We consider different geometries involving interfaces between both a monolayer and 
a bilayer and between two bilayers. The applied bias opens a sizable gap in the spectrum inside the 
bilayer barrier region, thus leading to large changes in the transmission probability and electronic 
conductance that are controlled by the applied bias. 
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I. INTRODUCTION. 

The idea of carbon based electronics has been around 
since the discovery of carbon nanotubes almost fifteen 
years ago. Much progress has been made but many 
problems associated with manufacturability remain still 
to be resolved - see for example the recent reviews in 
Refs. [HIH Recently, another possible platform for car- 
bon based electronics was discovered in graphene,— i.e., a 
two-dimensional (2D) honeycomb lattice of carbon atoms 
that can be viewed either as a single layer of graphite or 
an unrolled nanotube. The electric field effect has al- 
ready been demonstrated in these systems: by tuning a 
gate bias voltage one can control both the type (electrons 
or holes) and the number of carriers4^ For a recent re- 
view on the rise of graphene see Ref. [(J 

One fundamental difficulty with most of the graphene 
devices studied so far is the experimental fact that there 
exists a universal (sample independent) maximum in the 
resistivity of the order of 6.5kf2 near the Dirac point in 
all these systems. 4 - This relatively small resistivity lim- 
its the performance of devices via a poor on-off ratio. 
The reason for the nonzero minimal conductivity is the 
peculiar gapless spectrum and presence of disorder in 
the samples. Several theoretical studies using different 
methods find a universal minimum in the conductivity 
in some limits. But the reason for it and its exact value 
varies ^^ i 10 ! 11 ' 12 ^ 3 ' 14 see also the recent review Ref. [l5|. 
There has also been a number of earlier theoretical stud- 
ies of junctions and barriers in both monolayer and bi- 
layer graphene systemsi^ i 13 ' 16 i 17 ' 18 i 19 

One way of getting around the minimal conductivity 
is to use nano-ribbons made out of graphene. Because of 
the confinement these systems are generally found to be 
gapped, 20 which is consistent with recent experiments. 21 
In this paper we propose another simple geometry that 
transforms the semi-metallic graphene into a semiconduc- 
tor with a truly gapped spectrum without confinement. 
The presence of a gap makes the properties of the sys- 
tem more robust to perturbations. This is crucial for 
device performance and possible device integration since 
imperfections are always present. The basic idea is to 



use a bilayer region as a barrier for the electrons. By 
manipulating the electrostatics with gates and/or chem- 
ical doping a gap can appear in the spectrum. By tuning 
the chemical potential one can move the system from 
sitting inside of the gap, where the electronic transmis- 
sion is exponentially suppressed, into the allowed band 
regions where the transmission is close to one. Other in- 
teresting features of the proposed geometry are that the 
gap size can be tuned with gates, allowing for an exter- 
nal control of the electronic properties, and the bilayer 
can be integrated as one of the components of a pure 
graphene based device^ 2 - Recently barriers made out of 
double-gated bilayer graphene have been fabricated and 
characterized in Ref. |23| . and their results are consistent 
with a gate-tunable gap. 

The paper is organized as follows: In Section [II] we in- 
troduce the model of the system that we are considering 
and give explicit expressions for the wave functions. In 
Section HTIl we match wave functions considering different 
geometries. In particular we study a monolayer-bilayer 
interface, a bilayer-bilayer interface, and two bilayer bar- 
rier setups. In the barrier setups we use either monolayer 
graphene or bilayer graphene in the leads. In Section ITVl 
we present results for the conductance in the different 
barrier setups. A brief summary and the conclusions of 
the paper can be found in Section |Vj For completeness 
we also provide some details about the self-consistent de- 
termination of the gap in Appendix [Al In Appendix [Bl we 
provide some details of how we compute the transmission 
amplitudes. 



II. MODEL 

A schematic picture of the system is shown in Fig. [T] 
The basic structure involves a graphene sheet extending 
to the left and to the right of a region where there is a sec- 
ond graphene sheet sitting on top of the first, this region 
is the bilayer barrier (BB) region. The whole structure 
is assumed to sit on top of a dielectric spacer insulating 
the graphene from a back gate. Later on (Sec. IIII E|) we 
will also consider a system where the regions to the left 
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FIG. 1: [color online] Geometry of the biased graphene bilayer 
barrier with monolayer leads. 



and right of the barrier are made out of graphene bilay- 
ers. We will refer to the two systems as having monolayer 
and bilayer leads respectively. 

To the left (L) and to the right (R) of the barrier 
the low-energy effective Hamiltonian (near the K point 
of the Brillouin zone) has the form of the 2D Dirac 
Hamiltoniani^i 



V a /vF k x + ik % 
: x -ik y V a /v F 



(1) 



where k = (k x , k y ) is the 2D momentum measured rela- 
tive to the K point, a = L, R and vp « 10 6 m/s is the 
Fermi-Dirac velocity (we use units such that K = 1 = Vf 
from now on). The Dirac Hamiltonian acts on a spinor 
representing the wave functions on the two sublattices: 
tp = (ipAli V'bi) T - In the following, we mostly work with 
the case where Vl = Vr = so that the energy is mea- 
sured with respect to the Dirac point in the monolayer re- 
gions. The spectrum of the Hamiltonian in Eq. ([T]) is then 
simply E±(k) — ±k, (k 2 = k x + where the plus (mi- 
nus) sign is associated with electron (hole) states. The 
low-energ y eff ective bilayer Hamiltonian has the form [see 
e.g. Refs. HUH: 



Hbb — 



f Vi q x + iq y tj_ 

q x - iq y Vi 

. t± V 2 q x - iq y 

V q x + iq y Vi 



(2) 



Here the two-dimensional momentum is q = (q x , q y ) and 
the corresponding spinor is ijj = (ipAi, ipBi, i>A2, V'B2) T - 
t± f=a 0.35 eV is the hopping energy between nearest 
neighbor atoms in different planes (i.e., Al and A2). The 
monolayer is connected to plane 1 in the bilayer. Solving 
for the spectrum of the Hamiltonian in Eq. J2]) one finds 
four energy bands given by 

E±, s (q) = (Vi+V 2 )/2 



± 



2 



,1^4(^ + ^)^2+4^ (3) 



where s = ± and V = Vi — V%. Thus for the two bands 
closest to the Dirac point (i?±._) the spectrum is gapped 
and has an unusual "Mexican hat" dispersion as was 
pointed out in Rcf. [25j. There exists other examples in 
the literature of materials with a similar dispersion, of- 
ten called "camel-back" -dispersion instead of "Mexican 
hat" . This feature can arise - like in our case - in the 
k • p approximation when two bands that are close to 
each other in energy are allowed to hybridize. Exam- 
ples of materials where a camel-back has been proposed 
include Tellurium, 2 ! GaP, 2 ^ and GaAs^ 

The crucial property for the structure proposed here 
is that we allow for different voltages on the two layers 
in the BB region: V\ +- V%. This possibility have been 
noted before ) 25 ! 30 but here we are exploiting this feature. 
These potentials can be created by a uniform electric 
field through the bilayer that generates a charge imbal- 
ance between the two layers. In a transport measurement 
the graphene is connected to electron reservoirs and can 
hence become charged. When voltages are applied to the 
gates and the graphene structure charge is redistributed 
to minimize the total electrostatic energy. The problem 
is basically that of a capacitor. Upon studying the prob- 
lem the bilayer must be viewed as a single unit since the 
planes are connected by orbital overlap. For example, the 
induced charge imbalance between the layers will screen 
the applied electric field; t± also works against the ap- 
plied field since it tends to equalize the densities in the 
two planes. 

A simple estimate of these effects is provided by 
a self-consistent Hartree theory [see Appendix [5] and 
Refs. SHI | . 

For an isolated uncharged infinite bilayer 
we find that the net effect is to replace the applied volt- 
age difference V by a smaller effective Vmf- For some 
reasonable parameters [V < tj_ and interlayer distance 
d sa 3.35 - 3.6 A 3 ^ 3 -], Vmf is down by a factor of or- 
der three compared with the bare value. In particu- 
lar for an experimentally accessible voltage drop of 90 V 
over 300 nm we find an effective voltage difference of 
Vmf ~ 40meV. This value might be improved upon 
incorporating a dielectric (e.g., SiO^) and/or using a 
thinner dielectric spacer. Thus it is not unreasonable 
to have Vmf ~ 100 meV. This estimate was done be- 
fore we became aware of the measurements reported in 
Refs. and Hf| In the first of these references the gap 
is measured in ARPES to be as large as 200 meV, but 
the charge densities are also quite large in their case 
(n ~ 1 — 6 • 10 13 cm~ 2 ). In the second reference the 
maximum obtained gap is estimated to be ~ 100 meV. 
The largest possible value of the gap is estimated to be 
~ 300 meV and is limited by the dielectric breakdown of 
the Si0 2 . 

Experimentally the bias can be controlled by different 
methods. The conceptually simplest and most flexible 
method is to use a back gate and a top gate (like in 
a dual gate MOSFET geometry), preferably using split 
gates and therefore allowing for different gate voltages in 
the monolayer and bilayer regions. It is worth to mention 
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that local top gates have already been successfully fab- 
ricated on single-layer graphene samplesj 36 ' 37 ! 38 ' 39 The 
field is progressing rapidly, since for example only a year 
ago no top gate had been reported in graphene systems. 
Moreover, very recently a top gate was also fabricated 
on bilaycr graphene and characteristics consistent with a 
gate-tunable gap were reported^ 3 - Another possibility is 
to change the chemical environment by depositing donor 
or acceptor molecules on top of the structure. 34 ^ 5 These 
act like dopants in a semiconductor and allow for in- 
dependent control of the bias and the chemical poten- 
tial. Note however that this method always introduces 
impurities into the system with potentially important 
consequences.— 

On symmetry grounds a more general Hamiltonian 
than the one in Eq. © is certainly allowed as discussed 
in Ref. for the case of graphite. For example, the 
couplings 73 and 74 associated with electron hopping be- 
tween carbon atoms that are not nearest neighbors in 
different layers, familiar to the graphite literature, are 
possible. 32 In the BB the effects of these terms are less 
important than for the low-energy features in graphite 
since the gap is a robust feature in the low-energy spec- 
trum of the BB. The electrostatic response in different 
sublattices within each layer is also likely to be different, 
that is, one should use Vai 7^ Vbi etc. If the applied 
voltage difference is not too large compared with t± the 
implications of this effect are presumably small. The ap- 
plied field and the pressure from anything on the top of 
the structure will probably affect both the interlayer dis- 
tance and the interlayer coupling. The details of the band 
structure including all the effects above is a very compli- 
cated problem that has yet to be studied. Nevertheless, 
we believe that the Hamiltonian in Eq. correctly cap- 
tures the main features of the spectrum, including the 
important formation of the gap in the spectrum near the 
Dirac point. In view of the large uncertainty in the pa- 
rameters involved, it is meaningless to pursue a more 
complicated model at this point. The inclusion of the 
other parameters introduce no principal problems, but 
the analysis becomes more complicated. 



A. The eigenvectors 

The normalized energy E eigenvectors of H a in Eq. ([1]) 
can be written as: 



1 



E-V* 

V2(E-V a ) \k x -ik v 

when the particles are "on-shell" : 

(E - V a 



) 2 — + ky. 



(4) 



(5) 



A simple way to generate the eigenvectors of Hbb in 
Eq. ^ is to note that the columns of the Green's func- 
tion G° in the bilayer region (G° = \E — Hbb\ ) are 
proportional to the eigenvectors if the relation between 



the energy E and the momentum q are "on-shell" . It is 
straightforward albeit tedious to generalize this approach 
to the more general Hamiltonians discussed above. Us- 
ing this we extract two different energy E eigenvectors of 
Hbb as 



wbb,ai 



(UE-V 2 f-ql-q^\(E-V 1 )\ 

((E-V 2 f-q 2 x -qS(q x -iq y 
t 



\ 



_(E-V 2 )(E-Vi) 
t ± (E-V 1 ){q x +iq y ) ) 



(6) 



and 



«BB,A2 



/ t x {E-V 2 )(E-Vi) \ 
t±{E - V 2 ){q x - iq y ) 
[{E-V 1 f-ql-ql]{E-V, 



{[(E-V,) 2 



1y] {.Qx + 



(7) 



These just differ by their overall normalization. It is 
straightforward to check that these vectors are indeed 
eigenvectors of Hbb hi Eq. ^ by direct substitution. 
The "on-shell" condition in the bilayer region reads: 



2{ql + ql) = {E-V 1 f + {E-V 2 f 



± V [(E - V,y -(E- V 2 y] + At\(E - V 1 )(E - V 2 ). 



(8) 



III. DIFFERENT GEOMETRIES 



In this section we compute the transmission amplitudes 
for different edges and geometries. By convention the 
incident wave is taken to arrive from the left side of the 
barrier and is transmitted to the right. 



A. Zig-zag termination with monolayer leads 

With our conventions the zig-zag termination of the 
barrier corresponds to cutting the strip along the y- 
direction. For simplicity we consider the case that the 
width W of the structure is large enough so that the 
boundary conditions in the transverse direction are irrel- 
evant. It is then convenient to assume periodic bound- 
ary conditions and use translational invariance and fix 
k y = q y to be a good quantum number in addition to the 
energy E. This assumption can be relaxed ^ 19 ' 42 When 
the width becomes small enough that the quantization 
of the transverse momentum becomes important the sys- 
tem becomes similar to a semiconducting nanotube with 
a finite radius. As discussed by Brey and Fertig, it is 
a good approximation to use the continuum description 
of a graphene ribbon if it is wide enough and provided 
that the proper boundary conditions for the continuous 
model are employed.— For the zig-zag edges one can work 
in the single valley approximation and the correct bound- 
ary condition is to take the wave- function to vanish in one 
of the sublattices (A2 on the left and B2 on the right in 
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our case) at the BB boundaries. We choose the energy E 
and k y — \E\ sin(0) so that there are propagating states 
to the left of the junction, with k x = Ecos(<j)). The four 
solutions of Eq. |5]) for q x inside the bilayer we denote by 
±q x i and ±q x2 - 



B. Monolayer-bilayer interface 

Consider the step geometry where there is no right 
end of the bilayer. Then one must only keep states that 
propagate to the right or decay as one moves into the 
bilayer. Thus, to the left we take the wave function to 
be: 



i>L = lv L ,+e ikxX + rv Lt _e 



(9) 



where r is the reflection amplitude, and vl± are the 
spinors associated with the sublattices given in Eq. (TJ|. 
In the bilayer one has, 



V'bb = a 1+ v B B,i+e lq:clx + a 2 +v B B.2+e 



(10) 



where ai(2)± are scattering amplitudes and «bb,i(2)± the 
respective spinors computed from Eq. ((6|) or Eq. (0. 
Matching of the wave functions only involves their con- 
tinuity. Because the associated differential equation is of 
first order this is sufficient to insure current conservation. 
Explicitly the boundary conditions are: 





= 0)| A1 


= VW^ 


= °)U 


(11a) 




= °)Ibi 


= V'BB {% 


= 0)Ibi» 


(lib) 


i>BB{x 


= °)L 


= 0. 




(11c) 



From this we compute [for details see Appendix [B] the 
transmission probability T(E,<f>) = 1 — |r| 2 and the an- 
gular averaged transmission probability 



T{E) = / —TIE, 

J-tt/2 1" 



(12) 



Some representative results are shown in Fig. There 
are two features that are apparent in the figure: (a) There 
is a small asymmetry between the angles ±0. Hence, a 
current of electrons without valley polarization leads to a 
transmitted current with a finite valley polarization. This 
is not the case at other boundaries, like a potential step 
applied to a graphene monolayer or a graphene bilayer J£ 
The breaking of the symmetry between the two Dirac 
points arises from the lack of time reversal symmetry, 
as we consider a current carrying state, and the lack of 
inversion symmetry, induced by the zigzag bilayer edge or 
the bias potential in the bilayer (for a general discussion, 
see Ref. 144 for a particular discussion see Appendix |B|) . 
As a result, the barrier discussed here can be used as 
a device which creates a valley polarized current^ (b) 
There is a clear asymmetry between positive and negative 
energies. This can be understood by noting that the 



monolayer is coupled to layer 1 in the bilayer. When the 
energy is tuned to V\ the weight on sublattice Al goes 
to zero [c.f. Eq. Consequently the current in plane 
1 is zero at that energy, and hence no current can flow 
into the bilayer. It appears as though the current goes 
continuously to zero as the band edge is approached. At 
the other edge (E ~ V 2 ) of the band gap the current in 
plane 2 is zero [c.f. Eq. ©], but the current can now 
flow in through the other plane. This fact is responsible 
for the sharp edge in the transmission amplitude at the 
conduction band edge. 
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FIG. 2: [color online] Transmission amplitudes in the 
monolayer-bilayer step geometry. Left: Energy dependence 
of T for V2 — Vi — 40meV and different values of V%. 
Right: Angular dependence of T{E,<j>) for Vi = lOmeV and 
V2 = 50 meV and different values of the energy. 

It is important to choose the momenta for right-movers 
in Eq. ([9]) and Eq. (fT0|) such that their group velocity 
v g = dE / dq x > 0. Otherwise one may erroneously con- 
clude that T < for some values of the energy, this is 
sometimes referred to as the Klein paradox4£ It is also 
worth noting that the actual charge distribution near 
the edge is a complicated problem that involves a self- 
consistent solution of the Poisson equation and the band 
structure, beyond the scope of this study. This may lead 
to corrections to the simple wave function matching we 
use here. Furthermore, it is known that edges can lead 
to interesting effects in graphene systems such as edge 
states and self-dopingji 



C. Biased bilayer barrier with monolayer leads 

Consider the barrier geometry of Fig. [1] We assume 
the length L of the bilayer region to be large compared 
to the lattice spacing so that the continuum model is 
applicable. In this geometry one also needs the wave 
function to the right of the barrier: 



i/jr = tv Ri+ e 



ik x (x — L) 



(13) 



where t is now the transmission amplitude. Inside the 
barrier one generally needs all momentum components: 



V'bb = ai+^BBa+e 1 
+ ai„w B B,i-e~ 



a 2 +VBB,2+e Ha 

<X2-^BB,2-e' 



(14) 
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In this case, in addition to the boundary conditions in 
Eq. (fTTj) there are also those at the right edge: 



the barrier. This is related to the modulo 3 effect found 
in the spectrum of an armchair graphene nano-ribbon in 
the nearest neighbor tight-binding approximation. 43 





= l )\ai 


= $bb(x 


= L )\av 


(15a) 






= L )\m 




= L )\m> 


(15b) 


E. Bilayer-bilayer step and Biased bilayer barrier 


^BB [x 


= L )\B2 


= 0. 




(15c) 


with bilayer leads 



It is a simple task to match the boundary conditions 
and in this case one finds six equations for the six un- 
knowns: t, r, &i+, at—, a.2+ , and a2_. The results 
for T(E, L, <f) — \t\ 2 for some representative parame- 
ters are shown in Fig. [3] The oscillations in the trans- 
mission amplitudes are due to the possibility of hav- 
ing resonances inside of the BB region. For example, 
in the right panel of Fig. [3] the two energies are chosen 
such that they have the same q 2 computed from Eq. ([5]) 
for the propagating mode. In one case (E — 5meV), 
\E\ sin(</>) = q y <C q and the resonances are pronounced 
with the distance between consecutive maxima approxi- 
mately given by Tr/q x ~ n/q ~ 25 mi. In the other case 
(E = — 45meV), siii(0) = q y ~ q and most of the 
resonance phenomena averages out upon performing the 
angle average because of the larger variation of q x . 




FIG. 3: [color online] Transmission amplitudes in the BB with 
monolayer leads. Left: Energy dependence of T(E, L, (j>) for 
different angles and L = 50 nm. Right: Length dependence 
of T(E, L) for different energies. In both the figures Vi = 
and V2 = — 40meV. 



D. Armchair termination 

The calculations for the armchair geometry is similar 
albeit more involved since one cannot work in the sin- 
gle cone approximation^ 3 . In this case, one must instead 
mix the valleys to be able to fulfill the boundary condi- 
tions so that the wave functions on both of the lattice 
sites in layer 2 vanish at the boundaries. This doubles 
the size of the matrix problem that must be solved. For 
example, to compute the transmission amplitudes in the 
BB geometry one must solve a system of twelve equa- 
tions. The calculation is a straightforward extension of 
the case above and although the shapes of the curves are 
not exactly the same there are no new features except for 
a parity effect associated with the number of unit cells in 



It is important to note that it is not necessary to have 
the BB region defined by actually cutting the second 
graphene sheet. Another possibility is to use a bilayer 
throughout, and to use a local top gate to create a lo- 
cal gap and hence a barrier in the bilayer. As we shall 
see, the characteristics of this type of junction is better 
than the one with the monolayer leads. In particular the 
oscillations in the transmission amplitudes and conduc- 
tivities are much smaller because the matching between 
two bilayers is usually better than between a monolayer 
and a bilayer. 

We will consider the simple case that the bilayer in the 
L and R regions is essentially gapless, we also assume 
that \E\ < t± for the incoming wave as this will likely be 
the experimental situation. Since the gapless bilayer is 
a special case of the gapped bilayer with no gap we may 
use the formulas and spinors from Section III Al directl 
with V\ = V2 = 0. This case was considered in Ref. 
and allows for some further simplifications of the spinors, 
but that is not necessary here. 

The incoming wave is taken to be a traveling wave with 
absolute momentum given by k 2 = \E\(t± + \E\). As in 
the above we define k x — ±/ccos(</>) and k y = fcsin(</>). 
One should also take care to define the sign of k x for 
the incoming wave so that the wave is a rightmover. We 
will denote the corresponding spinors by VBo,on where 
a = ± goes with exp(±ik x x) and denotes right- and left- 
movers respectively. To be able to fulfill the boundary 
conditions one must also consider the decaying modes^ 
which will have an imaginary value of the momentum in 
the x-direction: k x = ik x . One can show that the correct 

value is n x — ^J\E\{t±_ — \E\) + k 2 . The corresponding 

spinors we write as ubcmq m an obvious notation. 

The calculation proceeds exactly as in the other cases 
when one has identified the particular incoming propa- 
gating mode in the bilayer to the left: 



lp L = lVBO, 



_e tk * x + rv B o -e lkxX 



r vbo, 



(16) 



Now we can consider a step geometry where the incoming 
wave from the unbiased bilayer propagates into a biased 
bilayer. In this case the spinor in Eq. (fH))) should be 
matched with the one in Eq. (|10[) at x — 0. More details 
for this case are provided in Appendix [Bj Some repre- 
sentative results are shown in Fig. [U and these should 
be contrasted with the case of a monolayer-bilayer step 
in Fig. [21 Note that the asymmetry between ±0 is also 
present in this case. For the case of a bilayer barrier with 
bilayer leads the spinor to the right is 



i'R = tv B o,+e 



ik x {x—lS) 



+ t' VB0,iZ 



(17) 
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FIG. 4: [color online] Transmission amplitudes in the unbi- 
ased bilayer-biased bilayer step geometry. Left: Energy de- 
pendence of T for V<2 — V\ = 40meV and different values of 
Vi. Right: Angular dependence of T(E,cf>) for V\ = lOmeV 
and V2 = 50 meV and different values of the energy. 



The wave functions in Eq. (fTBj) and (|17[) should be 
matched with the one in Eq. (jT5J) to the left and to the 
right. In this case the correct boundary conditions is to 
take all of the components of the 4-spinor wave function 
to be continuous at the two boundaries of the BB re- 
gion. The transmission amplitude is in this case given by 
T(E,(p) = \t\ 2 . An example of the transmission ampli- 
tudes is shown in Fig. [5] It is clear that the magnitude of 
the oscillations in the amplitudes are much smaller than 
in the cases involving monolayer leads. Moreover, the 
pronounced resonances found in the length-dependence 
of Fig. [3] is largely gone in this case. This is probably 
due to the fact that the wave functions of two bilayers 
are better matched than those of a monolayer and a bi- 
layer. It is also interesting to note that the effective gap 
becomes larger than the actual gap in the BB for larger 
values of the angles. This is due to the fact that, given 
the energy, the absolute value of the momentum is much 
larger in the bilayer than the monolayer. Consequently 
one has to go to larger values of the energy in Eq. ([8]) to 
have a mode that is not decaying inside of the BB region. 
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FIG. 5: [color online] Transmission amplitudes in the BB with 
bilayer leads. Left: Energy dependence of T(E, L, (f>) for dif- 
ferent angles and L = 50 nm. Right: Length dependence of 
T(E, L) for different energies. In both the figures Vi = and 
V 2 = -40meV. 



IV. RESULTS FOR THE CONDUCTANCE 

Using the Landauer formula [see e.g. Ref. |47[, we find 
that the current across the BB is given by 

/= 7T / dE"£\t n (E)\ 2 [f(E-» R )-f(E-n L )]. (18) 

n 

Here / is the Fermi distribution function and [i^ (/ir) is 
the chemical potential in the left (right) lead, n labels the 
modes and t n (E) the corresponding transmission ampli- 
tude at energy E. From this the finite temperature linear 
response conductance can be computed as a function of 
the overall chemical potential ji: 



G(/x) = — — / dEM(E)T{E) 



-af(E-v) 



dE 



(19) 



For the BB with monolayer leads, M(E) ~ W\E\/n is the 
number of transverse propagating modes in the mono- 
layer at energy E. For bilayer leads the relation is instead 
M(E) ~ Wy/\E\(t_i_ + \E\)/tt. At zero temperature the 
expression in Eq. (|19p simplifies to 



G(m) = ^M( M )T(/i). 



(20) 



Some results are presented in Fig. [5] and Fig. [7] As ex- 
pected, if the barrier is wide enough and the energy is 
tuned to be inside of the gap, the conductance is strongly 
suppressed. Outside of this region the conductance is 
larger by many orders of magnitude. It is also clear from 
Eq. (fT9|) that a finite temperature T will lead to a smear- 
ing of any sharp feature on an energy scale of approxi- 
mately AT. Nevertheless, as long as the temperature is 
much smaller than the gap large on-off ratios are possible. 

Let us finally comment on possible effects associated 
with roughness or impurities at the edges of the sample. 
These will induce some intervalley scattering and lead to 
an angle average. This may be a serious problem for the 
proposals which emphasize the angular dependence of the 
transmission?^ Because we are considering the transmis- 
sion integrated over the angles this should only weakly 
affect our results. The resonances are also likely to sur- 
vive in a real sample when many incoming modes overlap 
with one of the eigenmodes of the BB region. Roughness 
at the ends of the BB will probably broaden the reso- 
nances however. 



V. SUMMARY AND CONCLUSIONS 

We have studied the problem of electronic transmis- 
sion through a graphcnc bilayer barrier as a function of 
applied voltage between the layers and the overall chem- 
ical potential. We have considered two types of devices, 
one with monolayer leads and one with bilayer leads. In 
the first type the barrier region is defined by having a 
bilayer only in a small part of the sample. In the second 
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FIG. 6: [color online] Left (right): Semi- log (linear) plot of 
the zero-temperature conductance divided by the width W 
(in nm) of the BB device with monolayer leads in units of 
Go = 4e 2 /h as a function of the chemical potential [i. Vi — 
Vi = 40 meV and L = 50 nm 



0.16 
0.14 
0.12 
0.1 



C3 

g 0.08 
O 0.06 
0.04 
0.02 













Vi = - 20 meV 






Vi = -10meV 




















'/ : 








ft 









p (eV) 



1i (eV) 



FIG. 7: [color online] Left (right): Semi- log (linear) plot of 
the zero-temperature conductance divided by the width W (in 
nm) of the BB device with bilayer leads in units of Go = 4e 2 jh 
as a function of the chemical potential fi. V2 — Vi = 40 meV 
and L — 50 nm 



between the planes is not equal to the bare externally 
applied potential difference Vq. A simple approxima- 
tion that takes into account the screening of the external 
field by the BB is to include the interaction among the 
electrons within the BB at the Hartree level in a self- 
consistent manner. Such a calculation was applied for 
the half-filled case for the preprint of the present paper, 
and it was worked out independently at the same time for 
a more general case by McCann in Ref. Hfl More recent 
works along these lines includes a joint experimental- 
theory paper^ and ab-initio calculations. 31 There has 
also been a study of the screening of an external elec- 
tric field in graphene multilayer systems in the A-B 
stacking 4£ 

In this appendix we provide details of the self- 
consistent determination of the gap. In particular we give 
explicit analytic expressions for all of the quantities in- 
volved in the non-linear equation that needs to be solved. 
These expressions are potentially useful for anyone who 
wish to apply this simple theory. 



1. Hartree Approximation 



Charges on the two planes in the bilayer leads to an 
electrostatic energy given by the capacitive coupling 



E c = —2irde 2 Sn\n2, 



(Al) 



type of device the barrier is instead defined by a local 
gate in a system made entirely out of a graphene bilayer. 
The latter system seems to have a smoother electronic 
characteristics due to the absence of sharp boundaries 
that are present in the first device. We have shown that 
the transmission probability and the electronic conduc- 
tance are strongly dependent on the applied bias leading 
to the possibility of applying this geometry for carbon 
based electronics. 
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APPENDIX A: SELF-CONSISTENT 
DETERMINATION OF THE GAP 

It is important to note that due to the polarization 
of the BB the actual size of the potential difference V 



in the simplest approximation of two uniformly charged 
planes. First we decouple the contributions from the 
total density n — ri\ + n% and the density difference 
Sn = 112 — n\. Here we are interested in the latter term, 
which we will treat in the Hartree mean field approxi- 
mation. Denoting by (Sti(Vmf)} the expectation values 
of Sn in the ground state of the mean field Hamiltonian 
[i.e., Eq. © with the substitution Vi — Vi = V — > Vmf], 
the mean field equation can be written as 

V MF = V - 2nde 2 (5n(V M F)) , (A2) 

where the expressions for (Sn) at half-filling is given be- 
low in Eq. (|A3[) . If all the energies are expressed in 
eV the mean field equations becomes simply Vmf ~ 
Vb — 7.3 (<5ti(Vmf))- The solution of the mean field equa- 
tions at half-filling is shown in Fig. [5] Away from half- 
filling one must also include the corrections due to the 
partial filling of the conduction or valence band. 

An important quantity for the self-consistent determi- 
nation of the gap size is the density difference between 
the layers Sn. In the simple model we are using this quan- 
tity depends on t±, V and the density n. At half-filling, 
when the chemical potential is sitting inside of the gap 
the result to leading order in the cut-off (see below) is: 



8 



6n 



V 



2n(tl + V 2 ) 3 / 2 



VJt 2 + V 2 



(V 2 /2 + t 2 ± ) + V^tl + V 2 /4) 



+ t\{t\ + U 2 /2)ln 



t\ + V 2 /2 + y/tl + V 2 /A^t\ + V 2 

(yffiTV*-V)V/2 



(A3) 



The occupation asymmetry 5n as a function of V at half 
filling is also depicted in Fig. [8] One can easily check 



which include both spin polarizations and the two valleys. 
Now one may use that 




FIG. 8: Left: Occupation asymmetry 8n at half filling as a 
function of the bias V. Right: Self-consistently determined 
value of the bias potential Vmf as a function of the applied 
potential Vq. 



that the expression in Eq. (|A3|) reduces to the correct 
expression in the limit of decoupled planes (t ± = 0) . 



A^A2 + iVB2-^Ai-iVBi = -2V(uj 2 + k 2 -V 2 /4-ti/2), 

(A6) 



and 



D'(E ±!0 ) = (32E ±t ^ ±(V* +t\)k* +t\, (A7) 



to write 



Sn 



V 



d{k 2 ) 



x [E+,+ - fc 2 + F 2 /4 + t 2 x /2] . 



(A8) 



2. Occupation asymmetry at half-filling 

First we introduce the shorthand N a j/D = G° J QJ for 
the diagonal components of the bare Green's function 
that is defined by G° - r " - ^— 1-1 
Det[w- Wbb], and N, 
the matrix [u — 7Ybb]- Using this the density of states on 
sublattice aj can be written as 



= [u - Wbb]" 1 - We take D 
aj to be the appropriate cofactor of 



N, 



5(lu - E / 



(A4) 



Here the /3-sums are over the different bands and we have 
suppressed the frequency- and momentum-dependence of 
the functions for brevity. If we want to compute the 
density difference between the layers 8n = ri2 — n\ when 
the chemical potential is sitting inside the gap we must 
calculate 



Sn= did [p A 2 + PB2 - PA1 - PBl] 

J —oo 



A d 2 k Na2 + Nb2 — Nai — Nbi 
(2n) 2 



D'(uj) 



J=E- lf3 

(A5) 



Using E+.+ E+,- = y/(fc 2 - V/2) 2 + VH\ /4 one can 
convince oneself that the integral in Eq. (1A8|) is conver- 
gent as A — > oo so that the leading term is independent of 
the cutoff. Changing the integration variable to z defined 
by 



^4(V 2 +t 2 ± )k 2 +t\., 



(A9) 



the integral can be performed analytically with the result 
shown in Eq. (|A"3)) . 



APPENDIX B: EXPLICIT MATRIX EQUATIONS 



In this appendix we show how to obtain the matrix 
equations that we then solve numerically to obtain the 
transmission amplitudes. For the simplest case of a 
monolayer-bilayer step, using Eqs. ©, (fTU|) . and (JB]) 
the boundary conditions in Eq. |[TT|) can be rewritten as 
a matrix equation: 
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( E - \ 




h x iky 1 


■( 








[{E-V 2 f-ql l -kl]{E-V 1 ) [{E-V 2 f-ql 2 -., ll 

[(E-V 2 ) 2 -q 2 xl -k 2 ](q xl -ik y ) [{E-V2) 2 -ql 2 -kl]{q x 2-ik y )\ [ Oi+ ) . (Bl) 




For this simple case it is also possible to work out an explicit expression for r: 

_ _E 3 - 2V 2 E 2 - [hi + q xl + q 2 x2 - V 2 + q xl q x2 - k x (q xl + q x2 )] E — (k x — ik y ){q xl + q x2 )Vi 
E 3 - 2V 2 E 2 - [k 2 + q 2 xl + q 2 x2 - V 2 + q xl q x2 + k x (q xl + q x2 )\ E + (k x + ik y )(q xl + q x2 )Vi ' 

I 



(B2) 



Now one can substitute the correct values of the momenta 
[c.f. Eq. ©J such that the incoming state is a right-mover 
and that only states that decay or propagate to the right 
inside of the bilayer are present. We note that for energies 
such that q x \ is real and q x2 is imaginary - which is often 
the case for energies such that V < \E—(Vi+V 2 )/2\ <t± 
- the transformation k y — > —k y is in general not simply 
associated with a phase of r. Consequently there is an 
asymmetry in the transmission amplitude as shown in 

Fig. m 

The reason for the asymmetry is the broken inversion 
symmetry. Either the symmetry is broken by a zig-zag 
edge of the bilayer or by the bias field. It is clear that the 
bias potential breaks the inversion symmetry in the point 
in the middle between the two A atoms of the unit cell. 
This is crucial as it breaks the sublattice symmetry that 
is otherwise present. It is also important that there exists 
a mode that is evanescent as this breaks the symmetry 
between states with k y and — k y when one is matching the 



wave functions. For if all momenta are real, one can easily 
convince oneself that to solve for — k y one must only take 
the complex conjugate of the solution with k y , thus the 
only difference is the phase between r and r* which will 
not affect the transmission amplitude. In addition it is 
necessary that there exists a mode that can transmit the 
current. Otherwise all of the incoming current is reflected 
and no asymmetry can be generated, note that this is the 
case for a monolayer with a gap. 

The other more complicated cases involves 6x6 and 
8x8 matrices. The procedure to obtain the transmission 
amplitude is a straightforward generalization of the ex- 
ample worked out above, but the full form of the matrices 
are too long to write out here. The matrix inversion is 
then performed numerically. 

Another interesting example is to match a wave coming 
in from an unbiased bilayer and propagating into a biased 
bilayer. In this case the matrix equation can be written 
as: 



/ 




E [(E - V 2 f q 2 xl k 2 } (E V,) [{E V 2 ) 2 q 2 x2 - k 2 } (E - V,) \ 

-k x -ik y ) (~iK x -ik y ) [(E - V 2 f - q 2 xl - k 2 ](q xl - ik y ) [(E — V 2 ) 2 - q x2 — k 2 ] (q x2 — ik y ) 



E\ \E\ t±(E-Vx)(E-V 2 ) 

\sign(E)(-k x + iky) sign(E)(-iK x + ik y ) t±(E - Vi)(q x i + ik y ) 



t x (E-Vi)(E-V 2 ) 
t±{E-Vi)(q x2 + ik y ) J 



(B3) 



Here we have used the fact that the spinors simplify in 
the leads where the inversion symmetry is not broken (i.e. 
Vi =V 2 = 0).— Also in this case one finds that there is an 
asymmetry between negative and positive angles within 
each valley. In this case this is due to the fact that the 



inversion symmetry is broken by the bias potential in the 
BGB. In the case of Vi = V 2 ^ the inversion symmetry 
is not broken and the transmission is symmetric between 
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